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Abstract 

Background: Specialized interactions help structure communities, but persistence of specialized organisms is 
puzzling because a generalist can occupy more environments and partake in more beneficial interactions. The 
"Jack-of-all-trades is a master of none" hypothesis asserts that specialists persist because the fitness of a generalist 
utilizing a particular habitat is lower than that of a specialist adapted to that habitat. Yet, there are many reasons to 
expect that mutualists will generalize on partners. 

Plant-soil feedbacks help to structure plant and microbial communities, but how frequently are soil-based symbiotic 
mutualistic interactions sufficiently specialized to influence species distributions and community composition? To 
address this question, we quantified realized partner richness and phylogenetic breadth of four wild-grown native 
legumes {Lupinus bicolor, L arboreus, Acmispon strigosus and A. heermannii) and performed inoculation trials to test 
the ability of two hosts (L bicolor and A. strigosus) to nodulate (fundamental partner richness), benefit from 
(response specificity), and provide benefit to (effect specificity) 31 Brodyrhizobium genotypes. 

Results: In the wild, each Lupinus species hosted a broader genetic range of Bradyrhizobium than did either 
Acmispon species, suggesting Xhd^X Acmispon species are more specialized. In the greenhouse, however, L bicolor 
and A. strigosus did not differ in fundamental association specificity: all inoculated genotypes nodulated both hosts. 
Nevertheless, A. strigosus exhibited more specificity, i.e., greater variation in its response to, and effect on, 
Bradyrhizobium genotypes. Lupinus bicolor benefited from a broader range of genotypes but averaged less benefit 
from each. Both hosts obtained more fitness benefit from symbionts isolated from conspecific hosts; those 
symbionts in turn gained greater fitness benefit from hosts of the same species from which they were isolated. 

Conclusions: This study affirmed two important tenets of evolutionary theory. First, as predicted by the Jack-of-all- 
trades is a master of none hypothesis, specialist A. strigosus obtained greater benefit from its beneficial symbionts 
than did generalist L bicolor. Second, as predicted by coevolutionary theory, each test species performed better 
with partner genotypes isolated from conspecifics. Finally, positive fitness feedback between the tested hosts and 
symbionts suggests that positive plant-soil feedback could contribute to their patchy distributions in this system. 

Keywords: Mutualism, Symbiosis, Specialization, Coevolution, Lupinus bicolor (dove lupine), Acmispon strigosus 
(strigose trefoil), Lotus, Alpha-proteobacteria, Legumes, Rhizobia 



Background 

Specialized biotic interactions contribute to processes that 
structure communities [1] and adapt populations to local 
partners and environments [2,3]- But, why do specialists 
exist if a generalist organism can occupy more environ- 
ments and partake in more beneficial interactions [4,5]? 
The "Jack-of-all-trades is a master of none" hypothesis 
asserts that specialists persist in a heterogeneous 
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environment because the fitness of a generalist utilizing a 
particular habitat is lower than that of a specialist adapted 
to that habitat [5-7]. This hypothesis is based on the as- 
sumptions that adapting to new habitats involves fitness 
costs and that traits adaptive in one subset of environ- 
ments are negatively genetically correlated with traits 
adaptive in other environments [4,8-11], leading to local 
adaptation [2,3]. When the habitat is a host or interaction 
partner, this trade-off could lead to co-adaptation [12]. 

Fitness feedbacks between plants and specialized soil- 
based microbes significantly affect plant productivity. 



© 2014 Ehinger et a!.; licensee BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative 
Commons Attribution License (http://creativecommons.Org/licenses/by/2.0), which permits unrestricted use, distribution, and 
reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain 
Dedication waiver (http://creativecommons.0rg/publicdomain/zero/l.O/) applies to the data made available in this article, 
unless otherwise stated. 



Ehinger et al. BMC Ecology 2014, 14:8 
httpy/www.biomedcentral.com/l 472-6785/1 4/8 



Page 2 of 18 



community composition, and the distribution and abun- 
dance of plants [13-21]. Strong plant-soil feedbacks are 
often negative; due to antagonistic interactions between 
plants and pathogens [13,21-27]. However, many plants 
engage in mutualistic interactions with soil-borne symbi- 
onts [6]. Is there sufficient specialization among plant 
hosts to soil-borne mutualistic symbionts to impact plant 
community structure and diversity via positive feedbacks? 

Multiple arguments predict that mutualists interacting 
with partners acquired from the environment are likely to 
be generalists [28]. Since mutualists have higher fitness 
when partnered than when not, any partner should be bet- 
ter than no partner. If so, then selection should disfavor 
specialized mutualists (i.e., those with reduced interaction 
breadth relative to the breadth of available partners) that 
resist associating with common partners [29]. Similarly, 
specialized mutualistic lineages should be vulnerable to 
perturbations in partner populations [30-36] and the geo- 
graphic distribution and spread of specialized mutualists 
could be limited by partner availability [37-44]. General- 
ized mutualists might also access a wider range of envir- 
onmental conditions by utilizing partners with differing 
ecological tolerances [45-48]. Finally, evolutionary conver- 
gence of traits among mutualists within an interaction 
group [49-51], but see [52], could reduce variance in part- 
ner benefit and thereby weaken the effectiveness of selec- 
tion to specialize [53,54] . 

Specialized mutualists nonetheless exist, e.g. [55-57]. In- 
deed, meta-analysis of bipartite interactions found that mu- 
tualistic webs are more specialized than antagonistic ones 
[1]. Further, mutualism theory predicts that, when available 
partners vary in quality, selection favors mechanisms by 
which mutualists can select partners [58-63], which can 
lead to specialization. However, if more-beneficial partners 
are rare or distributed unpredictably, then selection for 
specialization via partner choice might be weak [45,64,65]. 
How frequently symbiotic mutualistic interactions are spe- 
cialized is therefore an open question. 

Here, we examined specialization in a group of wild le- 
gumes that interact with Brady rhizobium genotypes in 
coastal sand dunes of Sonoma County, California. Bradyr- 
hizobium is a genus of rhizobia that can infect legume 
roots and fix nitrogen endosymbiotically within root nod- 
ules in an interaction that is generally mutualistic [66,67]. 
Dominant genes that either restrict nodulation with or 
alter the effectiveness of particular rhizobial genotypes 
have been described in both natural and managed plant 
populations [57,68], but the effect of such genes on fitness 
feedbacks has been little studied. Moreover, little is known 
about whether legume- rhizobium interactions are suffi- 
ciently specialized to contribute to structuring communi- 
ties via plant-soil feedbacks or coevolution. 

We quantified the realized richness and phylogenetic 
breadth of rhizobial symbionts collected from plants 



growing in situ, and measured the frequency with which 
each host species associated with each of the identified 
Bradyrhizobium genotypes. Association frequencies quan- 
tify the interaction [69] or link [70] strength of each pair- 
wise interaction and were used to calculate the Paired 
Differences Index PDI [71], for each host. This index sum- 
marizes variation in a partner s link strengths [69,70] with 
all available partners without making any assumptions 
about their statistical distribution. 

Realized association specificity [72,73] is strongly in- 
fluenced by the local sampling environment [74]. There- 
fore, we experimentally paired partners under controlled 
greenhouse conditions to examine fundamental associ- 
ation frequency [4,75] of two host species. We also used 
the greenhouse experiment to quantify the fitness effect of 
each host on the inoculated Bradyrhizobium genotypes, 
which is a measure of the hosts potential impact on the 
rhizobium population and therefore its functional role 
within the community [76]. 

Finally, we used these data to test whether this system 
supports the "Jack-of-all-trades is a master of none" hy- 
pothesis, as well as its underlying assumption that, as fre- 
quently interacting partners become co- adapted, they 
become less adapted to other partners and consequently 
exchange weaker benefits. 

Results 

Isolate collection and sequence datasets 

We isolated Bradyrhizobium bacteria from nodules of 
four native California legumes growing on the Bodega 
Dunes at Bodega Marine Reserve and Sonoma Coast State 
Park, Sonoma County. Two species, Acmispon strigosus 
and Lupinus bicolon are fast-growing annuals that are 
patchily distributed across the Bodega Dunes. In contrast, 
A, heermannii is a decumbent suffrutescent perennial and 
L, arboreus is a large upright perennial shrub, both widely 
distributed across the dunes. 

Two DNA regions {NifD and ITS) were PCR ampli- 
fied and sequenced from two to eight isolates from 
each plant, which yielded 84 sequenced isolates, including 
81 NifD amplicons and 82 ITS amplicons (Genbank acces- 
sion numbers in Additional file 1: Table SI). Sequences 
that differed by one nucleotide were considered distinct, 
resulting in 22 ITS genotypes, of which 13 occurred in 
multiple nodules, and 22 NifD genotypes, of which 13 oc- 
curred in multiple nodules. There were 39 unique geno- 
type combinations. 

Realized in situ association specificity 
(realized niche breadth) 

The observed PDIs of the four hosts {L arboreus: 0.264, A, 
strigosus: 0.301, A, heermannii: 0.395, L bicolor: 0.433) all 
differed significantly from the expected joint PDI (0.069 ± 
95% confidence limit of 0.002) and from each other, which 
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indicates that each of the four hosts specialized on a subset 
of the available NifD types, but differed in degree of 
specialization; L bicolor was most specialized. 

Genetic breadth of symbiotic partners 

In both gene networks, Lupinus hosts associated with a 
genetically broader range of bacteria than did Acmispon 
hosts. Isolates from within a host genus tended to cluster 
together (Figure 1) and exhibited little overlap between 
host genera. Host genera shared no genotypes in the NifD 
network (Figure lA). A single common ITS genotype (SI) 
was shared between host genera and one genotype from 
Lupinus (S8) clustered closely to most genotypes from 
Acmispon (Figure IB). 

Genetic distance matrices among bacterial communities 

Bacterial populations differed significantly between the 
two host genera but not between species within a genus. 
For both DNA regions, significantly large Ogt values and 
corrected between-population mismatch rates occurred 
between Bradyrhizobium populations isolated from dif- 
ferent host genera, but not between populations isolated 
from different host species within a genus (Table 1). 
However, bacteria isolated from Lupinus were more gen- 
etically variable than those isolated from Acmispon. For 
both DNA regions, the average pairwise sequence diver- 
gences between Bradyrhizobium genotypes from Lupi- 
nus were two- to five-fold greater than they were 
between genotypes from Acmispon (Table 1). 

Distribution of bacterial genetic variance 

In a two-level hierarchical AMOVA (Table 2), neither 
ITS nor NifD sequence explained a significant compo- 
nent of variance between host species within genus. 
Most genetic variance occurred among plants within 
species (50.29% for ITS and 40.84% for Nip) and be- 
tween host genera (49.14% for ITS and 57.92% for 
NifD), Permutations suggest that the variance compo- 
nent due to host genus was not significant (Table 2); 
even though genera harbored distinct genotypes in both 
networks (Figure 1) and all inter-genus values of pairwise 
genetic differences and Ogt were statistically significant for 
both DNA regions (Table 1). This discrepancy probably 
arises from a lack of power in the AMOVA, which treats 
genus as a random effect, even though it is associated with 
only one degree of freedom. 

Fundamental association specificity 
(fundamental niche breadth) 

In the greenhouse inoculation experiment, neither test 
host exhibited fundamental association specificity; each 
could nodulate all tested Bradyrhizobium genotypes 
(Figures 2 and 3). 



G X G interaction (Response specificity) 

On average, inoculation increased host shoot dry weight 
over that of uninoculated control plants. The average 
increase was 2.8-fold for A, strigosus (Fi, sis = 64.261, 
p value < 0.0001) and 2.9-fold for L bicolor (Fi, 314 = 
387.191, p value < 0.0001). However, bacterial genotypes 
differed significantly in their effects on shoot dry weight 
(Figure 4, Additional file 2: Table S2 and Additional file 3: 
Table S3). Notably, several genotypes failed to improve 
host growth beyond that of uninoculated controls (17 for 
A. strigosus and three for L. bicolor, using the Tukey- 
Kramer honestiy significant difference test; 13 for A. stri- 
gosus and two for L, bicolor, using the less conservative 
Students t-test). 

The conservative analysis presented above does not test 
for between-host differences in response specificity, which 
is typically detected as a significant interaction between 
test host species and Bradyrhizobium genotype G x G 
interaction [77], in an ANOVA including data from both 
host species. In a less conservative analysis, host by geno- 
type (G X G) interactions were prominent in the best fit- 
ting model (model 16 in Additional file 4: Table S4). A 
significant G x G interaction of test host with Bradyrhizo- 
bium groups isolated from different host species (test host 
by species of origin interaction, ^3,25 = 19.03, p>F< 
0.0001; Additional file 5: Table S5), occurs because vari- 
ation in response to different Bradyrhizobium genotypes 
was large in A. strigosus and almost absent in L. bicolor 
(Additional file 6: Table S6). Further, the G x G interaction 
due to test host species by Bradyrhizobium genotype nested 
within species of origin (covariance parameter estimate ± 
1 S. E. = 0.07235 ± 0.024; Additional file 7: Table S7) 
accounted for 7.24% ± 2.28% of the variance in plant dry 
weight. This significant interaction indicates that even 
within a Bradyrhizobium group isolated from the same 
host species, genotypes differed in the symbiotic benefits 
they provided to the two different test host species. 

As expected from the conservative analysis, test host 
species differed significantly in dry weight (least square 
means [logits] of shoot dry weight ± 1 S. E. for A. strigo- 
sus = -4.80 ± 0.10, L bicolor = -2.40 ± 0.11; tg = -18.36, 
p < 0.0001) and groups of Bradyrhizobium genotypes iso- 
lated from different host species differed significantly in 
average effect on test host shoot dry weight (species of ori- 
gin effect ^3,25 = 24.56, p>F<OmOl; Additional file 8: 
Table S8). In particular, shoot dry weight produced by ge- 
notypes isolated from L bicolor differed from that produced 
by genotypes isolated from other species (Additional file 6: 
Table S6 and Additional file 8: Table S8). When averaged 
across the test hosts, there was no significant variance in 
shoot dry weight due to inoculation with different Bradyr- 
hizobium genotypes isolated from the same host species 
(strain(sp_orig) covariance parameter estimate ± 1 S. E. = 
0.01929 ± 0.020; Additional file 7: Table S7). 
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T4 



Figure 1 Genotype networks generated by the neighbor-net algorithm. (A) NifD locus, based on the distance calculated by MODELTEST, 
using a GTR + I model of sequence evolution [121]. (B) ITS region, based on the distance calculated by MODELTEST using a HKY+ I + G model of 
sequence evolution [122]. Dotted lines represent the position of a genotype on the network. Each box represents one genotype, each symbol 
represents a host species (filled circles represent bacteria harvested from nodules of A. strigosus; open circles, A. heermannii; filled triangles, L 
bicolon open triangles, L arboreus), and numbers indicate the number of strains representing a particular genotype; labels Tl to T31 in italics 
represent the Brodyrhizobium genotypes used in the greenhouse experiment (as listed in Additional file 1: Table SI). One genotype, T30, appears 
only in the ITS network, because we were unable to sequence its NIfD locus. Genotypes Tl 1 and T12 do not figure in either network because 
they were addenda from another study, used for reference, but not part of the original sample. 



The fold increase in growth of inoculated plants 
over uninoculated controls was strongly affected by 
the interaction of test host with source host (Figure 5; 
Xd.f. = 1 = 25.61, p < 0.0001). Further, the planned con- 
trast comparing the fold increase in fitness obtained 



from Bradyrhizobium genotypes isolated from conspe- 
cifics versus heterospecifics was highly significant (likeli- 
hood ratio }^ f = 2 = 28.76, p <5 x 10"^). Both main effects 
were statistically significant (test host species Xd.f. = i = 
5.059, p < 0.03; host species of origin xif. = i = 5.059, p < 
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Table 1 Analyses of population structure of 
Bradyrhizobium sampled from in situ nodules 



(a) ITS region 





A. strigosus 


A. ineermonnii 


L orboreus 


L bicolor 


A. strigosus 


6.6601 


0.0422"' 


0.5814* 


0.4405* 


A. heermannii 


0.4506"' 


13.1696 


0.5459* 


0.4126* 


L arboreus 


24.6912* 


23.4825* 


26.2134 


0.0058"' 


L bicolor 


18.1281* 


17.9398* 


0.0078"' 


35.922 


(b) nifD locus 




A. strigosus 


A. ineermonnii 


L orboreus 


L bicolor 


A. strigosus 


5.1895 


-0.0040"' 


0.5649* 


0.6024* 


A. heermannii 


-0.0217"' 


7.4971 


0.5830* 


0.6223* 


L orboreus 


21.6199* 


22.0285* 


25.1368 


0.0406"' 


L bicolor 


19.3541* 


19.6343* 


1.0162"' 


1 7.0395 



Significance, obtained from 1023 permutations: ns, not significant; *,p< 0.0001; 
significant values in boldface. 

Sequence divergence (along diagonal), corrected between-population differ- 
ences (below diagonal), and between-population Ost (above diagonal). 



0.03), but the mean growth effect on each host species 
changed rank between conspecific and heterospecific pair- 
ings (xif. = 1. = 2^,62, p = 4 X 10"^). 

Response specificity, as indicated by the paired differ- 
ence indices calculated from within-block growth re- 
sponses stimulated by the different Bradyrhizobium 
genotypes, differed significantly between host species and 
was more than three-fold greater in A, strigosus (mean 
PDI ± 1 S. E. = 6.6 ± 0.7) than in L bicolor (mean PDI ± 
1 S. E. = 1.7 ± 0.1). 



significantly among Bradyrhizobium genotypes on each 
host (Additional file 2: Table S2, Figures 2 and 3). 

On A, strigosus, over 30% of variance in total nodule area 
was due to Bradyrhizobium genotype: genotypes isolated 
from congeners produced significantly greater total nodule 
area than did those isolated from lupines. Ten percent of 
variance in nodule number was due to Bradyrhizobium 
genotype: genotypes isolated from congeners produced 
slightly fewer nodules than did those isolated from lupines. 
Seventy percent of variance in average individual nodule 
area was due to Bradyrhizobium genotype: genotypes iso- 
lated from congeners produced significantly larger nodules. 

On L bicolor, only 15% of variance in total nodule area 
was explained by Bradyrhizobium genotype. Although geno- 
types isolated from congeners produced significantly greater 
nodule area, the difference was slight. Bradyrhizobium geno- 
type explained almost 15% of variance in nodule number, 
but there was no significant difference between genotypes 
isolated from congeners versus Acmispon species. Less than 
20% of variance in the average individual nodule area on L 
bicolor was due to Bradyrhizobium genotype. Genotypes iso- 
lated from congeners produce marginally significantly larger 
nodules than did those isolated from Acmispon hosts. 

Paired difference indices {PDI) calculated from host ef- 
fect on Bradyrhizobium fitness under non-competitive 
conditions of single inoculations is presented in Table 3. 
When measured in terms of nodule number. A, strigosus 
was significantly more specialized than L, bicolor, whereas 
PDI of the two hosts did not differ when calculated in 
terms of total nodule area. 



G X G interaction (Effect specificity) 

Bradyrhizobium fitness was estimated from nodule 
number, total nodule area, and area of the average nod- 
ule. Both nodule number and total nodule area differed 



Association between bacterial genotypes and phenotypes 

The matrices of pairwise phenotypic differences in A, stri- 
gosus shoot dry weight and average nodule area were each 
significantly associated with pairwise genetic differences at 
the NijD locus; association with the ITS region was not 



Table 2 AMOVA of (a) ITS and (b) NifD sequences from Bradyrhizobium from in situ nodules 



(a) ITS region 



Sources of variation 


df 


Sum of squares 


variance component 


% of variance 


P 


Between genera 


1 


440.898 


10.5300 


49.14 


ns 


Between species w/in genus 


2 


26.594 


0.1229 


0.57 


ns 


Within species 


78 


840.667 


10.7778 


50.29 




Total 


81 


1308.159 


21.4308 






(b) NifD locus 


Sources of variation 


df 


Sum of squares 


variance component 


% of variance 


P 


Between genera 


1 


290.41 


10.0195 


57.92 


ns 


Between species w/in genus 


2 


20.031 


0.2149 


1.24 


ns 


Within species 


77 


416.246 


7.06457 


40.84 




Total 


80 


980.938 


1 7.2989 







Significance: ns, not significant; p < 0.001. 

Sampled host species: A. strigosus, A. heermannii, L bicolor and L arboreus. 
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Figure 2 Average number of nodules produced by each Bradyrhizobium genotype on each host species. Untransformed data. A) A. 
strigosus. B) L bicolor. Error bars represent ± 1 standard error. Solid grey bars represent bacteria originally isolated from A strigosus and A. 
heermannii nodules {A bacterial group); solid black bars represent bacteria harvested from L. bicolor and L orboreus nodules (/. bacterial group). 



significant for shoot dry weight and marginal for average 
nodule area. Traits measured on L bicolor exhibited no sig- 
nificant associations (Table 4). 

Fitness feedback between host and symbiont 

For A, strigosus, there was an overall positive relationship be- 
tween total nodule area and host dry weight (genotype mean 



correlation between total nodule area and host dry weight = 
0.77, prob(Ho: r = 0) < 0.0001). Bradyrhizobium genotypes 
that were more beneficial for A, strigosus also produced on 
that host fewer nodules (genotype mean correlation between 
nodule number and host dry weight = -0.35, p < 0.0001) that 
were larger (genotype mean correlation between average 
nodule area and host dry weight = 0.71, p < 0.0001). 
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Figure 3 Average nodule size produced by eacli Bradyrhizobium genotype on eachi hiost species. Untransformed data (cm^). A) A. 
strigosus. B) L bicolor. Error bars represent ± 1 standard error. Solid grey bars represent bacteria originally isolated from A. strigosus and A. 
heermonnii nodules {A bacterial group); solid black bars represent bacteria harvested from L. bicolor and L. orboreus nodules {L bacterial group). 



A weak positive relationship between total nodule 
area and host dry weight (genotype mean correlation 
between natural log of nodule area and host dry 
weight = 0.36, prob(Ho: r = 0) < 0.05) suggests a weaker 
but still positive fitness feedback between L. bicolor 
and Bradyrhizobium symbionts. In L. bicolor, there 
was no relationship between host dry weight and nod- 
ule number (genotype mean correlation between 



natural log of nodule number and host dry weight = 
0.15, prob(Ho: r = 0) > 0.4). 

Missing samples and cross-contamination 

Nodule data were lost from 14 L. bicolor plants. Seven 
out of 100 control plants had nodules, indicating cross- 
contamination (perhaps water draining from containers 
was splashed from the greenhouse bench). These plants 
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Figure 4 Mean shoot dry weight of each host species when inoculated by each Bradyrhizobium genotype. Untransformed data (g). 
(A) A. strigosus. (B) L bicolor. Error bars represent ± 1 standard error. Horizontal reference lines represent mean dry weight of uninoculated 
control plants. Light grey bars represent bacterial group A, originally isolated from A strigosus and A. heermannii nodules; dark grey bars represent 
bacterial group L, originally isolated from L bicolor and L arboreus nodules; solid white bars represent uninoculated control plants. Genotype 
order recapitulates that found in Figures 2 and 3. 



(two A. strigosus and five L bicolor) were removed from 
analyses. Cross-contamination probably occurred late in 
the experiment because the nodulated control plants 
were still very small and had formed only a few small 
nodules. We are confident that these potential cross 
contaminations do not undermine our conclusions, 
since any broader cross contamination would only re- 
duce phenotypic differences observed between Bradyrhi- 
zobium genotypes. 

Discussion 

Association richness and phylogenetic niche breadth 

When growing wild, all four legume species that we exam- 
ined specialized on subsets of the nodulating community 
of Bradyrhizobium genotypes. However, accounting for 
genetic breadth of symbionts in wild-collected nodules 



shifted the ranking of host specialization. When 
Bradyrhizobium genotype was ignored (i.e., by simply 
counting the total proportion of symbionts that were non- 
identical at the focal loci), L bicolor ranked as most 
specialized and L arboreus as most generalized. However, 
Bradyrhizobium isolates from wild-collected nodules were 
genetically clustered by host genus and each Lupinus 
species hosted a phylogenetically broader group of symbi- 
onts than did either Acmispon species. Clustering was 
strongest when isolates were categorized by NijD geno- 
types, which suggests that hosts respond to trait loci located 
on the Bradyrhizobium symbiosis island. The observation 
that species in Acmispon exhibit greater realized association 
specificity than do those in Lupinus agrees with previous 
surveys of wild-collected nodules from plants in these 
genera [66,78-80]. However, neither L bicolor nor A, 
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A. strigosus L. bicolor 

Source species 



Figure 5 Mean shoot dry weight response to Bradyrhizobium 
genotypes isolated from conspecific versus allospecific hosts. 

Response = fold increase over uninoculated control plants. Closed 
circles indicate growth of A. strigosus test plants; open circles indicate 
growth of L 6/co/or test plants. Error bars represent ± 1 standard error. 



strigosus exhibited fundamental association specificity in 
greenhouse inoculation tests: both hosts could nodulate all 
tested Bradyrhizobium strains. 

A gap between fundamental and reaUzed association 
frequency is not uncommon [6] and indicates that envir- 
onmental context strongly influences nodulation. The 
complex biotic community within natural soil could influ- 
ence a genotypes nodulation success [81-83]. For example, 
soil pathogens could trigger systemic acquired resistance 
in the host, which might affect nodulation by some rhizo- 
bial genotypes but not others [84]. 

The availability in soil of multiple rhizobial genotypes 
also provides scope for plants to actively prefer certain ge- 
notypes. Host-symbiont recognition signaling is a well- 
studied feature of legume-rhizobium interactions [85]. In 
particular, it has been hypothesized that, during nodule 
formation and prior to nitrogen fixation, plants can 
recognize and prefer more-beneficial genotypes via pre- 
infection partner choice; reviewed in [62]. It is unclear 
however, whether this host ability would be evolutionarily 
durable, as it requires genetic linkage between at least 
three sets of genes in two different organisms (nitrogen 
fixation genes and signaling genes in the bacterial genome 



and receptors to the bacterial signals in the host plant gen- 
ome). Such linkage could be stable under very specific 
conditions, but could also be broken by any incident of re- 
combination [86]. In several legume-rhizobium systems, 
including A, strigosus [87], multiple -genotype inoculation 
experiments have found that more-beneficial and/or co- 
evolved genotypes are more likely to occupy nodules 
[81,88,89]. However, none of these studies could defini- 
tively distinguish plant choice from interactions among 
rhizobia, which could be due to varying abilities to com- 
pete for soil resources [90], withstand direct interference 
by other genotypes [91], colonize the rhizosphere [92], 
and utilize a range of rhizosphere resources [93,94]. In ex- 
periments with Bradyrhizobium strains collected from our 
site, neither in vitro growth rate nor interstrain interfer- 
ence was correlated with nodulation rate on A, strigosus 
[87], but these and other hypotheses should be further 
tested in more complex conditions where bacterial traits 
not directly related to symbiont effectiveness might influ- 
ence nodulation ability. 

Realized association frequency might also be restricted 
by patchily distributed plant and symbiont genotypes 
[95-97], which could be caused by coevolution [98] and/or 
soil heterogeneity. In the Bodega Dunes, two episodes of 
Holocene dune advance [99,100] left contrasting soils that 
might structure both plant and Bradyrhizobium popula- 
tions: L bicolor is restricted to mid-Holocene dunes whereas 
A, strigosus and A, heermannii occur only on poorly stabi- 
lized late-Holocene dunes (E. L. S. and T. J. M., personal ob- 
servation). In contrast, L arboreus occurs across dunes of 
both ages [101]. We specifically sampled nodules from L 
arboreus across its habitat range to control for potential 
confounding of host and Bradyrhizobium distribution; in- 
deed, L arboreus nodules harbored genotypes from across 
each gene network. However, the genotypic composition of 
isolates from L bicolor did not differ signiflcantly from those 
isolated from L arboreus^ which suggests that soil habitat is 
not the only determinant of genotype in wild-collected nod- 
ules of Lupinus, 

Response and effect specificity, fitness trade-offs, and 
fitness feedback 

The "Jack-of-all-trades is a master of none" hypothesis was 
supported by evidence that specialist A, strigosus obtained 



Table 3 PDI values for two Bradyrhizobium fitness components 

Fitness L bicolor A. strigosus Difference 

component ^^^^ p^^, ^CL UCL Mean PDI LCL UCL Diff PDI LCL UCL 

Nodule area* 0.217 0.189 0.246 0.228 0.186 0.271 0.102 -0.034 0.056 

Nodule number^ 0.672 0.605 0.739 1.067 0.832 1.302 0.395 0.183 0.607 

*Natural log transformed. 
tSquare-root transformed. 

For each statistic, LCL = lower 95% confidence limit and UCL = upper 95% confidence limit; Diff PDI = PDI difference between species. Variables that differed 
significantly between hosts are in bold face. 
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Table 4 Mantel tests of the correlation between genetic 



distance matrices and phenotypic difference matrices 



r 


P 


(a) A. strigosus dry weight 


NifD 0.3190 




ITS 0.0659 


ns 


(b) A. strigosus average nodule area 


NifD 0.413 




ITS 0.128 




(c) L bicolor dry weight 


NifD -0.006 


ns 


ITS -0.038 


ns 


(d) L bicolor average nodule area 


NifD -0.0996 


ns 


ITS -0.0280 


ns 



Significance, estimated from 50,000 permutations: ns, not significant; 
p < 0.05; **, p < 0.01 ; p < 0.001 . 



more fitness benefit from its genetically narrower group of 
beneficial symbionts than generalist L bicolor obtained 
from its genetically broader group of compatible symbionts. 

As predicted by the co-adaptation hypothesis, each host 
received greater fitness benefit from genotypes isolated 
from congeners or conspecifics (genotypes isolated from L 
bicolor provided particularly poor benefit to A, strigosus). 
Although the sample sizes used to estimate fitness are not 
large, this pattern of response specificity has been detected 
in several studies of unmanaged legumes [102-104] and 
suggests that symbiont effectiveness may be increased by 
host-symbiont coevolution [45,89,98,105]. 

A, strigosus plants were highly specialized in their fitness 
effect on the tested Bradyrhizobium genotypes, which is 
consistent with previous evidence that this host imposes 
absolute sanctions [64] on less effective genotypes [87]. 
There was much less variance among Bradyrhizobium ge- 
notypes in the benefits they obtained from L bicolor. 

Also supporting the co-adaptation hypothesis, genotypes 
of Bradyrhizobium were better adapted to hosts related to 
those from which their ancestors were isolated. This pat- 
tern was strongest on A, strigosus: genotypes isolated from 
congeneric hosts produced greater total nodule area and 
therefore likely yielded more progeny than did genotypes 
isolated from lupines. 

As a consequence, symbiont benefits to A, strigosus fit- 
ness positively fed back to symbiont fitness via increased al- 
location to nodules: genotypes beneficial to this host 
realized greater fitness benefits than did those that were 
not. An overall positive correlation between NifD genotypic 
mean values of total nodule area and host shoot dry weight 
on A, strigosus suggests that strong positive fitness feedback 
between mutualists in this partnership [67] is driven by 
traits encoded in the symbiosis island. This feedback 



pattern was much weaker on the generalist, L, bicolor, and 
was not significantly correlated with ITS genotype. 

Positive partner feedback can produce positive fre- 
quency dependence and drive the most beneficial partners 
to local fixation [49,106], a process that has been predicted 
for other legumes that exhibit symbiont specificity [83]. 

Much stronger positive fitness feedback might allow 
specialist A, strigosus to outcompete generalist L, bicolor 
in local patches [107]. However, if L, bicolor were a better 
colonizer, which is expected of generalists [108], then it 
might persist in a spatially patchy environment through a 
competition-colonization trade-off [109]. In that case. A, 
strigosus and L, bicolor would not coexist at a local scale, 
but both might persist at a larger spatial scale via patch 
dynamics. This process could also spatially structure Bra- 
dyrhizobium populations [12]. Indeed, wild A, strigosus 
plants were found in association with a genetically narrow 
range of Bradyrhizobium genotypes that were most benefi- 
cial to their growth, suggesting that these genotypes dom- 
inate Bradyrhizobium populations in areas inhabited by 
this host. In contrast, as predicted from the weaker corre- 
lations between plant and rhizobial fitness components on 
L, bicolor, wild L, bicolor plants associated with a genetic- 
ally wider range of Bradyrhizobium genotypes that were 
not necessarily the most beneficial. 

Conclusions 

A field survey detected differences in partner specificity 
among four sampled host species and greenhouse experi- 
ments revealed that at this site A, strigosus clearly special- 
izes on a genetically narrower range of symbionts than does 
L, bicolor. The fitness benefits that these two hosts received 
from symbiotic partners affirmed two important tenets of 
evolutionary theory. First, as predicted by the Jack-of-all- 
trades is a master of none hypothesis, specialist A, strigosus 
obtained greater benefit from its genetically narrow 
group of beneficial symbionts than generalist L, bicolor 
obtained from a genetically more diverse group. Second, 
as predicted by the co-adaptation hypothesis [2,3,50], 
each test species performed better with partner geno- 
types isolated from conspecifics. Further, host fitness benefit 
translated directly into symbiont fitness via increased alloca- 
tion to nodules: nodules occupied by genotypes isolated 
from conspecific hosts received the greatest benefit. These 
patterns were strongest in A, strigosus. Positive fitness feed- 
back between this specialized host and its co-adapted symbi- 
onts could drive positive plant-soil feedback and contribute 
to natural distribution patterns observed in the field. 

Methods 

Collection of wild Bradyrhizobium isolates 

During the end of March and the beginning of April 2007, 
we isolated Bradyrhizobium bacteria from nodules of four 
native California legumes growing on the Bodega Dunes at 
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Bodega Marine Reserve (38°19'01"N, 123°04'18"W) and 
Sonoma Coast State Park (38°20'34"N, 123°03'32"W). Por- 
tions of these dunes were stabilized in the mid-twentieth 
century by planting European beach grass, Ammophila are- 
naria [110,111], but no samples were collected in areas oc- 
cupied by A. arenaria. The four host legumes belong to two 
genera, Acmispon (A heermannii Greene [112], and A. stri- 
gosus Nutt. Ex Torr. & Gray [112], previously Lotus heer- 
mannii and Lo. strigosus, respectively) and Lupinus {L 
bicolor and L arboreus), which are from distantly related 
clades of papilionoid legumes [113]. 

We harvested and cultured isolates from up to 20 nod- 
ules from each of four to six seedlings of each host species 
using the procedures of Sachs and colleagues [66]. Briefly, 
intact plants were excavated, the root systems washed in 
tap water; nodules excised, surface sterflized, rinsed, 
crushed and individually streaked onto two replica 
plates of solid modified arabinose-gluconate medium 
(MAG). Successful cultures were archived in 25% 
glycerol-MAG at -80°C. 

Molecular methods and analysis 

Genomic DNA was purified from 25 \A of each frozen iso- 
late using Gentra Puregene Yeast/Bacteria kits, (Qiagen, 
Valencia, CA). We PGR amplified the intragenic spacer 
between 16 s and 23 s ribosomal subunits {ITS^ 1256 nt) 
[114] and portions of the nitrogenase a-subunit gene 
{NifD, 756 nt) [115] as previously described [66]. Amplifi- 
cation products were sequenced in both directions using 
an Applied Biosystems 96 capillary 3730x1 DNA Analyzer 
(Foster City, CA) at the University of California, Berkeley, 
Sequencing Facility. Two to eight isolates from each plant 
were successfully amplified and sequenced, which ultim- 
ately yielded a sample of 84 sequenced isolates (Additional 
file 1: Table SI). 

Sequences were aligned using MAFFT [116] with de- 
fault parameters. Genotypes were identified in MacClade 
4.05 [117]. Some portions of the ITS included indels and 
these regions were removed from the analysis. 

Quantifying realized in situ association specificity 

Specialization is most easily quantified as the number of 
taxa with which a taxon interacts [118], termed "partner 
richness". However, this measure accounts for neither as- 
sociation frequency [72,73], phylogenetic distance among 
partners, e.g., phylogenetic breadth [28], nor fitness effects 
of interactions. Resurging interest has stimulated new 
methods for quantifying ecological specialization [1,71,75]; 
we used the Paired Differences Index. 

From field coUected samples, we estimated the in situ 
link strength of each pairwise interaction, Pip as the pro- 
portion of nodules sampled from the host species oc- 
cupied by the Bradyrhizobium NifD type: 



(1) 



where n^ = the number of nodules on the i^ host species 
occupied by the Bradyrhizobium NifD type and ni = the 
total number of nodules sampled from the i^ host (viache = 
19, nAcsT= 18, niuAR = 21, nLUBi= 23), 

Observed specificity of each of the four hosts was then 
measured as the differential frequency of interaction 
with the 22 different avaflable NifD types, using the 
Paired Differences Index (PDI) measure [71]: 



PDIi 



(2) 



where is the strength of the strongest link with the 
i^^ host. The PDI values were not normally distributed; 
therefore, to compare PDIs of the four hosts, we esti- 
mated confidence limits from 1000 simulated 20-nodule 
populations generated as follows. Using the RAND 
("TABLE",) function of SAS® 9.3, we specified the fre- 
quency with which each genotype appeared in each simu- 
lated population, assuming that hosts randomly sampled 
rhizobia and therefore the observed frequencies of the 22 
NifD genotypes in the total sample of 81 nodules repre- 
sented their availability in the Bradyrhizobium population. 
From each simulated population, we calculated the joint 
PDI for all four hosts together: 



PDI^ 



joint 



N-l 



(3) 



where Pj is the link strength averaged across the four 
hosts of Bradyrhizobium strain Pi is the average link 
strength of the Bradyrhizobium strain with the strongest 
link averaged over the four host species, and N=20 nod- 
ules, then used these values to estimate the mean and 
variance of PDI and calculate 95% confidence limits. 

Genetic breadth of partners isolated from 
wild-grown hosts 

We constructed a molecular network for each DNA re- 
gion with the neighbor-net algorithm implemented in 
SplitsTree 4.8 [119], using the Akaike Information Criter- 
ion (AIC) implemented in Modeltest (version 3.4) to 
choose a suitable model of sequence evolution [120]. The 
NifD genotype network was estimated using the GTR + 1 
model [121] of sequence evolution and the ITS network 
was estimated with the HKY + I + G model [122]. When 
there were multiple isolates with the same sequence, we 
randomly picked one for analysis. 

A molecular network approach was favored over a phylo- 
genetic tree approach, as previous studies have repeatedly 
found incongruence between phylogenies of different gene 
loci in Bradyrhizobium bacterial strains [66,79,123-125]. 
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Moreover, in our sample, SplitsTree found significant re- 
combination in the ITS region and marginally significant 
recombination within the NifD locus. Genetic networks are 
standard tools for the analysis of populations in which ge- 
nomes recombine; such events will induce a reticulated 
representation of genealogies, which violates the bifurcating 
model of phylogenetic trees, as reviewed in [126]. 

For each DNA region, we computed the matrix of 
average number of differences between pairs of se- 
quences (pairwise differences) within and between host 
genera to test whether bacterial communities harbored 
by the two genera differed in genetic diversity. We 
further computed matrices of within and between host- 
species pairwise Ogt [127] to test whether pairs of host 
species harbored genetically distinct groups of bacteria. 
The statistic for nucleotide diversity is analogous to 
Wrights FsT) but accounts for sequence divergence as 
well as genotype frequency. Finally, each DNA region 
was analyzed with a two-level hierarchical AMOVA 
model with the following levels: within host species, 
between species within host genus, and between host 
genera. Significance tests for each level of the AMOVA, 
the values, and the pairwise differences were ob- 
tained by permuting genotypes between levels or popu- 
lations. These analyses were performed with Arlequin 
3.0b [128]. 

Greenhouse experiment 
Seed source 

Seeds of A. strigosus and L bicolor were obtained from 
ripe fruits of many plants growing at Bodega Marine Re- 
serve in May-June 2008, pooled, mixed and stored dry at 
room temperature. In mid-July, seeds were surface- 
sterilized and scarified [66,129], then germinated individu- 
ally in wells of 96-well culture plates containing 200 ml of 
sterile ddH20 (incubated at 15°C, 5-7 days). Seedlings of 
A, strigosus were transferred to bleach-sterilized 38-mm x 
140-mm Conetainers (Stuewe & Sons, Corvallis, OR, 
USA) filled with autoclaved quartzite sand [66]. One-week 
old seedlings of L, bicolor were transferred to bleach- 
sterilized 64-mm x 250-mm Deepots (Stuewe & Sons, 
Corvallis, OR, USA) filled with autoclaved calcine-clay 
(Turface MVP^ Profile Products, Buffalo Grove, IL, USA). 
Media were chosen to maximize plant survival. Neither 
medium provides nitrogen; we have no reason to expect 
that differences in media would affect symbiont specificity 
or fitness benefits to hosts. 

All transplants were placed in a greenhouse under ca. 
50% shade for seven days of hardening, with misting twice 
daily, and subsequently exposed to full sunlight, watered 
daily, and fertilized weekly with 10-ml Jensenss nitrogen- 
free solution [130] until treatments were applied one 
month after germination (see below). 



Preparation of bacterial inocula 

For the greenhouse inoculation experiment, we chose 29 
Bradyrhizobium strains that represented the range of a 
preliminary concatenated nifD and ITS network of the 
sequenced isolates. For reference to previously described 
Bradyrhizobium genotypes, we added two strains (Til 
and T12) from among those described by Sachs and col- 
leagues [66], for a total of 31 strains. A culture of each 
strain was initiated from -50 \A of original stock archived 
in glycerol, plated onto solid MAG medium [78] and incu- 
bated at room temperature until bacteria covered the 
plate. Bacteria were washed from plates into 50-ml poly- 
propylene tubes containing sterile dd H2O and vortexed. 
Bacterial concentration of each inoculum was estimated 
via light absorbance at 600 nm with a spectrophotometer 
and a previously established standard growth curve; then 
adjusted to 10^ cells mP^ by dilution in sterile dd H2O. 
Each treatment plant received 10 ml of bacterial suspen- 
sion and control plants received 10 ml of sterile dd H2O. 

Experimental design 

The two host species were planted and inoculated at the 
same time, but grown in different substrates and separated 
on the greenhouse bench to prevent shading of A, strigo- 
sus by the larger L, nanus. Each host species was arranged 
into ten 36-seedling blocks: 31 test seedlings, each ran- 
domly assigned to a Bradyrhizobium strain, and five seed- 
lings that received the sterile water control inoculum. All 
treatments were spatially randomized within each block 
and the 10 blocks were randomly located on the green- 
house bench, for a total of 360 plants per host species. 
Seven to eight weeks after inoculation, plants were har- 
vested to obtain oven-dried aboveground (shoot) dry 
weight (a component of plant fitness). We counted, ex- 
cised, and photographed nodules with a stage micrometer 
to obtain their total projected area using ImageJ64 
[131]. Because symbiotic rhizobia reproduce clonally 
and are ultimately released from legume nodules into 
the soil, nodule number and mass or projected area 
have been treated as multiplicative components of rhi- 
zobial fitness [67,129,132-134]. When obtained using a 
standard-density, single -genotype inoculum, nodule num- 
ber indexes the probability that an individual rhizobium 
cell can initiate a nodule [88,135]. Nodule size, measured 
as biomass or projected area, is a good proxy for the num- 
ber of viable rhizobia in a nodule [129,132,136]. Since 
most nodules are initiated by a single rhizobium cell 
[88,137], nodule size indexes rhizobium fecundity. Like 
many fitness components, nodule number and individual 
nodule size are often negatively correlated [138]. Nonethe- 
less, when all nodules on a host are occupied by a single 
genotype, the ratio of total nodule area to nodule number 
estimates average fecundity on that host for cells of that 
genotype. 
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Analysis 

Fundamental association specificity 

No analysis was required because all hosts nodulated all 
Bradyrhizobium genotypes. 

Host response specificity 

Because the experimental design technically precludes treat- 
ing test host as an independent variable in the ANOVA, it 
was most conservative to analyze each host species separ- 
ately. These models (Additional file 2: Table S2) included 
random effects of spatial block and Bradyrhizobium geno- 
type, with genotype nested within host species from which 
it was isolated and further nested within host genus of iso- 
lation, along with all appropriate interactions. A planned 
one degree-of-freedom contrast compared the growth of 
inoculated versus uninoculated test plants and a posterior 
1-df contrast compared the growth of plants inoculated 
with bacterial genotypes that had been isolated from the 
each of the two plant genera (/.e., Lupinus, L, versus 
Acmispon, A), Bacterial genotype effects on dry weight 
(Additional file 3: Table S3) were compared using a 
Tukey-Kramer test of honestly significant pairwise dif- 
ferences [139]. 

Although the experimental design technically pre- 
cludes treating test host as an independent variable in 
the ANOVA, if the different growth conditions (differ- 
ent growth media and locations in the greenhouse) of 
the two test host species did not interact with the effect 
of Bradyrhizobium genotype, then both species might 
be analyzed with a single model. Lack of replication of 
Bradyrhizobium genotype within blocks precluded test- 
ing for this interaction within each species. Nonetheless, 
we cautiously tried this less conservative approach, in 
which a significant interaction between test host species 
and Bradyrhizobium genotype G x G interaction [77], 
would indicate differences in response specificity be- 
tween L. bicolor and A. strigosus. 

We analyzed the entire data set (both host species) by 
adding the fDced factor of test host and all its interactions to 
a nested hierarchy of generalized linear mixed models to 
examine how host species, Bradyrhizobium genotypes, the 
species or genus of plant from which they were isolated, and 
two-way interactions affected plant benefit. The two hosts 
differed so much in response specificity that they exhibited 
very different distributions of shoot weight. For this reason, 
all models employed a gamma probability density function 
and a logit link function. Corrected Akaike Information Cri- 
teria (AICc) were used to find the best model (i.e. the one 
with the lowest AICc). Relative likelihoods of other models 
were compared using ^(^^^^^-^^1)12^ ^^^^^ AlC^in is the 
corrected AICc of the best-fitting model and AICi is the 
AICc of the model being compared [140]. This analysis was 
performed with SAS Version 9.3 (SAS Institute Inc., 2002- 
2010, Gary, NC, USA). The best fitting model (model 16 in 



Additional file 4: Table S4) included the fDced effects of host 
species (host), host species from which the Bradyrhizobium 
genotype was isolated (sp_orig), and the host species by spe- 
cies of origin interaction, random effects of block (block) 
and Bradyrhizobium genotype nested within the host spe- 
cies from which the genotype was isolated (strain(ge_orig)), 
and the interaction of host species with Bradyrhizobium 
genotype nested within host species of origin (host* strain 
(sp_orig)). Thus, the G x G interaction was split between a 
fDced effect interaction (Additional file 5: Table S5) and a 
random effect interaction (Additional file 7: Table S7). 

We also compared test host growth response, relative 
to uninoculated control, to Bradyrhizobium genotypes 
isolated from conspecifics versus response to those iso- 
lated from allospecifics. To calculate growth response 
(fold increase over uninoculated control), we used shoot 
weight of the control plants in each block to calculate the 
mean growth effect on each test host of the genotype in 
each block: {Wi - Wo)/wo, where Wi = shoot dry weight of 
the plant inoculated with genotype / and Wq = shoot dry 
weight of the uninoculated conspecific control plant. This 
variable was analyzed with a generalized linear model in- 
cluding test host, host species of origin, and their inter- 
action, with an exponential distribution and canonical 
(reciprocal) link function. A planned contrast was used to 
compare the fitness benefit obtained from Bradyrhizobium 
genotypes isolated from conspecific versus heterospecific 
hosts. This analysis was performed with JMP Version 
10.0.0 (SAS Institute Inc, Gary, NC, 1989-2007). 

To quantify specialization using the host fitness re- 
sponses, we calculated the fitness effect of each Bradyr- 
hizobium genotype on each test host species separately 
for each of the k = 10 blocks. 



Pijk 



^ijk-^iCk 



(4) 



where Wij^ = shoot dry weight of the ith host species in- 
oculated with the yth strain in the /<th block, and Wick = 
shoot dry weight of the uninoculated conspecific control 
plant in the same block. The Paired Differences Index 
(PDI) [71] of the ith test host can then be calculated for 
each block by ranking within-block values of pip^ from 
highest (j = I) to lowest (j = M), where N=30 tested 
strains, to obtain: 



PDIn 



Ef=2{PW<-Pijk) 



N-1 



(5) 



Variance among blocks was used to calculate 95% con- 
fidence intervals for the PDI of each test host species. 

Specificity of host effect on symbionts 

The effect of the two test hosts on the fitness of each bac- 
terial genotype was examined by univariate ANOVAs on 
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nodule number, total nodule area per plant, and average 
nodule size. All dependent variables were log transformed 
to meet the ANOVA normality assumptions [141], pp. 
185, 202-204. Standard deviations were not equal among 
factor levels; so we performed a Welch ANOVA, which al- 
lows for heteroscedasticity among factors [141], pg. 183. 

Specificities of effect of the two test hosts on Bradyr- 
hizobium genotype fitness under conditions of individ- 
ual inoculation were calculated from equation (4), 
above, by defining pijj^ in terms of each of two rhizobial 
fitness components, nodule number and total projected 
area of nodules. To meet the assumptions of normality 
and homoscedasticity between populations, nodule 
number was natural log transformed and nodule area 
was square root transformed. Variance among blocks 
was used to calculate 95% confidence intervals. 

Test of association between bacterial genotypes 
and phenotypes 

The groups of Bradyrhizobium genotypes isolated from 
congeneric versus heterogeneric hosts were identified 
by the population genetic analysis. The A group (geno- 
types Tl - T12) were originally isolated from Acmispon; 
the L group (genotypes T13 - T31) were originally iso- 
lated from Lupinus, Having been classified by host asso- 
ciation, the bacteria within one group may not be more 
related to each other than to bacteria in the other 
group. Thus, although this contrast compares the aver- 
age effects of the two bacterial groups, it does not test 
for a specific association between bacterial genotype 
and phenotype. To explicitly test for an association be- 
tween bacterial genotype and symbiotic phenotypes, we 
used Mantel tests of non-random association between 
matrices. We computed six matrices: two matrices of 
pairwise genetic distances between all pairs of geno- 
types, one for NifD and one for ITS, and four matrices 
of pairwise mean phenotypic differences between all 
pairs of bacterial genotypes, two for A. strigosus and 
two for L bicolor. Each genetic distance matrix was 
then compared with the four different phenotype matri- 
ces. Under the null hypothesis of these tests, a random 
association between the elements of the phenotype 
matrix and those of the genotype matrix would indicate 
that across Bradyrhizobium genotype pairs, the two ge- 
notypes within a pair shared phenotypes independently 
of the genetic distance between them. In contrast, a 
positive, non-random association would indicate that 
the phenotypic difference between the genotypes in a 
pair increases in concert with an increase in their gen- 
etic distance. The two reference genotypes. Til and 
T12, were omitted during computation of the matrices, 
as was genotype T30, which we were unable to sequence 
at the NifD locus. The ANOVAs were performed with 
JMP^ Version 7. (SAS Institute Inc., Gary NC, 1989 - 



2007); the Mantel tests were performed with the R stat- 
istical software [142]. 

Fitness feedback between host and symbiont 

To estimate fitness feedback between host and symbiont, 
the JMP statistical package was used to calculate the Pear- 
son product mean correlation between host fitness, mea- 
sured as shoot dry mass, and symbiont fitness, measured 
as nodule number, total nodule area per plant, and average 
nodule size. 

Availability of supporting data 

The sequence datasets supporting the results of this article 
are available from Genbank (http://www.ncbi.nlm.nih.gov/ 
genbank/); accession numbers JQ230720 - JQ230882. 

Additional files 



Additional file 1: Table SI. Title of data: Field-sampled nodules, their 
source host, identifyir^g symbols, and accession numbers for NifD and ITS 
sequences submitted to GenBank. Description of data: * host = host species, 
pin = plant within species, nod = nodule on plant, strain = strain identifier, 
T-type = inoculation type in greenhouse experiment (blank cells indicate 
strains not used in that experiment), g'type NifD and g'Xype ITS identify the 
NifD and ITS genotypes, respectively, included on Figure 1. Respective NifD 
and ITS GenBank sequence IDs and accession number are indicated as seq 
nifD, seq ITS, accession # NifD, and accession # ITS. Blank cells in these 
columns indicate isolates that did not amplify. 

Additional file 2: Table S2. Title of data: Separate ANOVA for each host 
species of plant and bacterial fitness-components. Description of data: 
Significance: ns, not significant; **, p<0.01 ; p<0.001 . ^ Degrees of 
freedom for F-test of random block effect: numerator = 9; denominator for 
(a) = 271, for (b) = 256. ^ Degrees of freedom for F-test of random genotype 
effect: numerator = 30, denominator for (a) = 271, for (b) = 256. ^ Tests of 
contrast between Bradyrhizobium groups A (isolated from Acmispon) and /_ 
(isolated from Lupinus); degrees of freedom: numerator = 1, denominator for 
(a) = 271, for (b) = 256. 

Additional file 3: Table S3. Title of c/ata.- Tu key-Kramer honestly 
significant differences in trait values between Bradyrhizobium 
genotypes inoculated on two host species, (a) A. strigosus, (b) L. bicolor 
Description of data: Genotypes sharing the same capital letter are not 
statistically significantly different at p < 0.05. * shoot dry weight. 

Additional file 4: Table S4. Title of data: Generalized linear mixed 
models fitted to variance in host shoot dry weight. Description of data: ^ 
estimated G-matrix not positive definite. Host = test host species, block = 
spatial block, strain = Bradyrhizobium genotype, sp_orig = host species from 
which Bradyrhizobium genotypes were isolated, ge_orig = host genus from 
which Bradyrhizobium genotypes were isolated. All models assumed a 
gamma distribution with a logit link function. Model 16 (bold face) provided 
the best fit (i.e., minimized AlCc). Adding the block by host species 
interaction improved the fit (e(^i^nnin-Aici)/2 ^ ^ x 10-13) but the proportion 
of variance explained by this component (4.91% + 2.42%, Additional file 7: 
Table S7) was only marginally significant. The next best fitting model (model 
1) included the fixed effect of host and the random effects of block and 
Bradyrhizobium genotype, with interactions of host with strain and with 
block. Residuals plots suggested that this model did not perform much 
worse than the best model, but nesting Bradyrhizobium genotype into host 
species from which it was isolated greatly improved the relative likelihood 
that the model fit the data (model 16 versus model 1 ; e^^'^^'^^^'^'^^^ = ^ 28 x 
10-12)). Examination of Q-Q plots revealed that all other models provided 
much poorer fits to the data. 
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Additional file 5: Table S5. Title of data: ANOVA table of Type III tests 
of fixed effects on test host shoot dry weight. Description of data: From 
model 16 of Additior^al file 4: Table 54. 

Additional file 6: Table S6. Title of data: Differences of least square 
means (logits) of shoot dry weight of test hosts by host of origin. 
Description of data: From model 16 in Additional file 4: Table 54. Test host 
species: AC5T = A.strigosus; LUBI = L bicolor. Origin host species: ache = 
A. heermannii] acst = A.strigosus; luar = /_. arboreus; lubi = L bicolor. 

Additional file 7: Table S7. Title of data: Covariance parameter 
estimates for random factors associated with test host shoot dry weight. 
Description of data: From model 16 of Additional file 4: Table 54. 

Additional file 8: Table S8. Title of data: Least squares means (logits) of 
shoot dry weight of test plants inoculated with genotypes isolated from 
the four hosts. Description of data: From model 16 of Table 3. 
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